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ABSTRACT 

We use a set of large, high-resolution cosmological iV-body simulations to examine 
the redshift-space distortions of galaxy clustering on scales of order 10 — 200ft.~^Mpc. 
Galaxy redshift surveys currently in progress will, on completion, allow us to measure 
the quadrupole distortion in the 2-point correlation function, ^(cr, tt), or its Fourier 
transform, the power spectrum, P{k, /i), to a high degree of accuracy. On these scales 
we typically find a positive quadrupole, as expected for coherent infall onto overdense 
regions and outflow from underdense regions, but the distortion is substantially weaker 
than that predicted by pure linear theory. We assess two models that may be regarded 
as refinements to linear theory, the Zel'dovich approximation and a dispersion model 
in which the non-linear velocities generated by the formation of virialized groups 
and clusters are treated as random perturbations to the velocities predicted by linear 
theory. We find that neither provides an adequate physical description of the clustering 
pattern. If used to model redshift space distortions on scales for 10 < A < 200/i~^Mpc 
the estimated value of /9 (/? = f{no)/b where /(r^o) ~ ^^o'^ and b is the galaxy bias 
parameter) is liable to systematic errors of order ten per cent or more. We discuss 
how such systematics can be avoided by i) development of a more complete model of 
redshift distortions and ii) the direct use of galaxy catalogues generated from non- 
linear A^-body simulations. 

Key words: cosmology: theory - large-scale structure of Universe - galaxies: clus- 
tering - galaxies: distances and redshifts. 



. ^ ' 1 INTRODUCTION voids and other underdense regions. Thus the anisotropy of 

' , , , , . ^ , , , ^ galaxy clustering encodes information about the galaxy ve- 

, in early galaxy redsnitt surveys (eg. Jackson 1972: Gregory , r. u , n j i • j c u ^i, ^ 

5h „ j a j j \ a loj locity tield and hence the underlying mass density held that 

rri ' & Thompson 1978), some or the most striking artefacts ■ j_ im j_ r i i j- j_ j_- 

'■^ ' , ,„ r ^ 1, ? gave rise to it. Measurements oi redshiit-space distortions 

observed were the so-called nngers oi God , ridges m the r i i . i , • , .i i •. 

can thereiore be used to place constraints on the density 

galaxy distribution pointing directly at the observer. What j_ r-. i .li i • ^ 7 i • i i ^ 

" 10 .J ri jT_-r parameter, ilo, and the bias parameter, 0, which relates the 



galaxy distribution to the mass distribution. 



had been observed was the first evidence of the redshift- 
space distortion of galaxy clustering. If galaxy motions were 

perfectly described by pure Hubble flow, then a galaxy's red- Interesting constraints on the combination /3 = n^'^/b 

shift would be an accurate indicator of its distance and the have been obtained by analysis of optical redshift surveys 

clustering pattern observed in galaxy redshift surveys would (Loveday et al. 1996; Ratcliffe et al. 1997a) and IRAS sur- 

be statistically isotropic. The observed anisotropic clustering ^gys (eg. Hamilton 1993; Fisher et al. 1994; Cole, Fisher & 

pattern is generated by galaxy peculiar velocities which per- Weinberg 1995 and more recent papers such as Fisher & 

turb galaxy redshifts and hence their inferred distances. On Nusser 1996; Bromley, Warren & Zurek 1997). However, the 

small scales, the non-linear velocities of virialized groups and statistical errors remain large and no consensus on the value 

clusters create the 'fingers of God' by stretching out these of (3 has been reached (for a review see Strauss & Willick 

structures along the line of sight. On large scales the cluster- 1996). The next generation of large surveys, including the 2- 

ing pattern is predicted (as suggested by Sargent & Turner degree-Field (2dF, CoUess 1996) and Sloan Digital Sky Sur- 

1977) to be compressed along the line of sight by coherent in- vey (SDSS, Gunn & Weinberg 1995), will enable redshift- 

fall onto galaxy clusters and super-clusters and outfiow from gpace distortions to be measured with unprecedented ac- 
curacy. These surveys should allow (3 to be measured to 
an accuracy of a few per cent, and possibly for Qq and b 

* S.J.Hatton@durham.ac.uk to be separately constrained. To achieve this accuracy and 

t Shaun.Cole@durham.ac.uk avoid systematic errors it is important that the theoretical 
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modelling of redshift-space distortions be accurate. This is a 
non-trivial requirement, as these surveys will provide their 
most precise measurements of the distortions on scales of 
30-100/i.~^ Mpc, where the predictions of pure linear theory 
are not expected to hold. 

In this paper we use a set of high-resolution Af-body 
simulations which accurately follow the evolution of cluster- 
ing over the range of scales which will be probed by the next- 
generation surveys. Rather than constructing mock galaxy 
catalogues wc instead use the full simulations to quantify 
the redshift-space distortions as accurately as possible, and 
through this assess the accuracy of the linear theory model 
and two proposed extensions to it. In Section 2 we define our 
notation convention and introduce the quadrupolc statistic 
which we use to assess the level of redshift-space distortion. 
Section 3 details the parameters of our set of A/^-body sim- 
ulations, and explains the methods we employ to construct 
biased galaxy samples. The linear theory prediction and the 
two other models are described in Section 4. In Section 5 
we compare the analytic models with measurements of the 
redshift-space distortions in the A'^-body simulations. We 
discuss results in Section 6. 



2 ANISOTROPY IN REDSHIFT-SPACE 
2.1 The power spectrum 

Our method of analysis uses the power spectrum, i-'(k), the 
Fourier transform of the correlation function. 



P(k) 



l^We^-d^r, 



(2.1) 



where we use the Fourier convention that the wavenumber, 
k, is the angular frequency corresponding to wavelength A, 
ie. k = 27r/A. The correlation function is given by 



e(r) = (5(x + r)5(x)), 

where 5{r) is the fractional overdensity, 



S{r) = 



P(r) - P 



(2.2) 



(2.3) 



The power spectrum has been used by many authors to 
study the growth of gravitational clustering in simulations, 
and to quantify the clustering observed in redshift surveys 
(for recent examples, see da Costa et al. 1994; Tadros & 
Efstathiou 1995; Lin et al. 1996). 

When galaxy distances are measured in redshift space, 
their peculiar velocities (ie. relative to the pure Hubble flow) 
distort the pattern of galaxy clustering by displacing galaxy 
positions along the line of sight. Thus a galaxy whose true 
position is r appears in redshift-space at the position 



• + U{r)i 



(2.4) 



where the line-of-sight peculiar velocity U{r) — v(r).f. Here 
we have adopted units in which the Hubble constant has unit 
value, we measure distance in units of h'^ Mpc, where h = 
i!fo/(100kms~'^Mpc~'^), and velocities in units of lOOkms"^. 

In redshift space this displacement in a preferred direc- 
tion causes the observed power spectrum to be anisotropic, 
with different values for wave- vectors along the line of sight 
to those perpendicular to it. Thus, the redshift-space power 
spectrum can be thought of as a function of k (= |k|) and /u, 



the cosine of the angle between wave-vector k and the line 
of sight, 



P(k) = P(fc,M). 



(2.5) 



In what follows wc adopt the convention that, by P(k) or 
P{k, n) we refer to the redshift-space power spectrum, since 
this quantity depends on both the magnitude and direction 
of k, whereas P(k) represents the real space power spectrum, 
depending only on the scalar k. 



2.2 The quadrupole ratio 

The redshift-space distortions can be conveniently quantified 
by a simple statistic. The anisotropy in P(k) is symmetric in 
/X, ie. P{k, ji) = P{k, —fi), so the distortion depends only on 
even powers of /i. To measure the extent of deformation from 
isotropy, we decompose the power spectrum into multipole 
moments using the Legendre polynomials, Li{n). Thus, 



P(fe,/i) = ^P(fc)L,(/x), 



(2.6) 



1=0 



where the sum is over even values of I. From the orthogo- 
nality relation for Legendre Polynomials, the multipole mo- 
ments Pi{k) can be found by evaluating 



Piik) 



P{k,n)Li{ix)^ 



(2.7) 



We find that estimates of Pi{k) rapidly become noisy for 
multipoles with Z > 2, so we choose to use as the key 
statistic for our analysis the quadrupole to monopole ratio, 
P2{k)/Po{k). 

We have chosen to carry out all our analysis in fc-space 

using power spectrum multipoles Piik). However, our results 
can readily be translated into predictions of the multipole 
moments of the correlation function, (,i{r) using the iden- 
tities derived in appendix B of Cole, Fisher & Weinberg 
(1994). 



3 SIMULATIONS 

In order to study the effects of non-linear gravitational evo- 
lution on redshift space clustering we examine a selection 
of iV-body simulations with different values of the density 
parameter f2o and the galaxy bias. The simulations we use 
are taken from a set of large volume simulations which have 
been produced in order to make mock redshift catalogues of 
the Sloan Digital Sky Survey (Gunn & Weinberg 1995) and 
the Anglo- Australian 2-Degree Field Survey (CoUcss 1996), 
as will be described in Cole, Weinberg, Hatton & Frenk (in 
preparation). The full set of simulations consists of two se- 
ries. The first is normalized to the amplitude of CMB fiuctu- 
ations in the COBE-DMR data set (Smoot et al. 1992) and 
the second normalized to match present day observations of 
galaxy clusters and large scale structure. Both series cover 
a range in f2o and include both open models (with Ao = 0) , 
and fiat models (f2o + Ao = 1, where Ao is the cosmological 
constant in units of 3-ffo). 

The main cause of redshift-space distortion on small 
scales is the random velocities of galaxies in groups and clus- 
ters. Since our aim is to study deviations from linear theory 
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Figure 1. Power spectra of the real-space mass distributions of 
the Af-body simulations described in Section 3. In each case the 
thin line is the linear theory, T = 0.25 CDM spectrum (3.1) 
which was used to set up the initial conditions of the simulations, 
the thicker line is the evolved, non-linear power spectrum of the 
simulation measured at z = 0. 



on scales of 10-200/i~^ Mpc this important rion- linear effect 
most be included at a realistic level. Therefore we have se- 
lected simulations from the series that are normalized to pro- 
duce approximately the observed abundance of rich galaxy 
clusters. This set of simulations have the rms mass fluctu- 
ation in spheres of radius 8/i~^ Mpc set to era = CSSJIq "'^ 
(White, Efstathiou & Prenk 1993). The shape of the initial 
power spectrum is the same in each simulation and is given 
by the Bardeen et al. (1986) formula for the cold dark matter 
spectrum, 



P{k) oc 



k X [ln(l + 2.34g)/2.34g]2 



[1 -I- 3.89g + (16.1g)2 + (5.46g)3 -|- (6.71g)4]i/2 



.(3.1) 



Here q = k/T and we set F = 0.25, as suggested by ob- 
servations of large scale structure (eg. Feldman, Kaiser & 
Peacock 1994; Maddox, Efstathiou & Sutherland 1996). To 
illustrate our results we focus mainly on two models, the 
Qq = 1 model and a low density f2o = 0.3, Aq = 0.7 model. 

The initial particle positions and velocities were set 
by using the Zol'dovich (1970) approximation to perturb 
particles from an almost uniform "glass" distribution. This 
"glass" distribution was produced using the technique de- 
scribed by White (1994), and Baugh, Gaztaiiaga & Ef- 
stathiou (1995). The simulations were evolved using the 
AP'^M code of Couchman (1991) using 192^ particles in 
a periodic cube of side Lbox = 345.6/i~^ Mpc (comoving 
distance). The softening parameter of AP^M's triangular- 
shaped cloud force law was set to r? = 270/i~^ kpc, 15 
per cent of the grid spacing. This choice corresponds ap- 
proximately to a gravitational softening length c = r//3 = 
90h~^ kpc for a Plummer force law. Further details and some 
analysis of these simulations can be found in Eke, Cole & 
Prenk (1996) and in Cole et al. (in preparation). 



3.1 Biased Galaxy Catalogues 

Results from recent galaxy surveys suggest that galaxies do 
not perfectly trace the mass distribution but rather pick 
out the highly non-linear regions. For instance. Peacock & 
Dodds (1994) show there are in fact different bias factors for 
differently selected galaxy samples, eg. IRAS galaxies and 
radio galaxies. Many of these recent investigations favour a 
value of /3 ~ 0.5, so for a flat universe with zero cosmological 
constant a significant amount of biasing is needed to match 
observations. To this end we employ three different methods 
of biasing our Qu = 1.0 simulation, in order to find out 
how robust our parameter estimation techniques are to the 
precise bias prescription. 

In each method we build up a 'probability field' on a 
grid, based on the density field smoothed on a particular 
scale. We then randomly sample the particles using the val- 
ues of the probability field at their nearest grid points £is 
their Poisson mean probability of being selected as galaxies. 

The first two methods we employ are Lagrangian, in 
the sense that the selection probability is based on the mi- 
tial density field. The third is an Eulerian one, where the 
galaxy distribution is 'painted on' depending on the final, 
evolved density field. In each case the bias is chosen to ob- 
tain a af^^, the variance in galaxy counts in spheres of radius 
8/i~^ Mpc, consistent with the results of Maddox, Efstathiou 
& Sutherland (1996), who find af' = 0.96 for the APM 
Galaxy Survey. Thus, the bias factor, defined as af^^ /as, has 
value h = 1.81. All our methods are local schemes based on 
the value of the overdonsity at a given point. The methods 
are expected to produce a constant bias on large scales. 

Method 1 High Peaks Biasing. For this method we 
smooth the initial density field by applying a sharp cut- 
off in the power at a scale rs, set to be, approximately, the 
shortest resolvable wavelength in the box, Vs « 4/i~^ Mpc. 
We then divide by the resulting mass variance, fjs, to ob- 
tain a new variable, !/(r) = S{r)/as. We use the results of 
Bardeen et al. (1986) to calculate the mean number of peaks 
above a certain threshold, Vt, for a random field, Gaussian- 
smoothed on a scale rgai = 0.54/i~^ Mpc, ie. the number of 
galaxy-scale peaks expected as a function of u (cf. White et 
al. 1987). The value of the threshold is adjusted until the 
required level of bias is obtained, from erf*' = bas- We find 
that Vt = 0.99 gives the correct normalization. 

Method 2 For this method we obtain z^(r) by smoothing 
the initial density with a Gaussian of scale rs — 5h~^ Mpc. 
Our biasing probability field is then a continuous function 
of u and two adjustable parameters, a and /3, given by 



P(i>) = A exp(a!/ -I- l3v\i>\) 



(3.2) 



where the normalization constant, A, is determined by the 

constraint 



/ 



2w 



(3.3) 



In the limit of small {5 <^ 1) large-scale perturbations in the 
initial density field, it can be shown that the bias on large 
scales is given by 



a 23 /"„,,,, exp( — ;y 12) , 
6=1-1 h— / P{v) \v\ ' ' dy. 



(3.4) 



■ ' ■ ■ ^/27r 

To fix a and (3 we demand that the rms galaxy density fluc- 
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tuations in two different sized cubic cells (5 and 20h^^ Mpc) 
match those predicted by assuming a power spectrum shape 
from the result for the APM survey found by Baugh & 
Gaztafiaga (1996). 

Method 3 In this model the final density field is 
smoothed on a scale rs = 5h~^ Mpc. The selection prob- 
ability field has a value or 1 according to whether the 
overdensity is below or above a threshold, t. We find that 
T = 1.53 provides the required bias. The particle distribu- 
tion resulting from this scheme is expected to contain large 
voids since the sharp cut-off at the threshold is extremely ef- 
fective at evacuating underdense regions (Weinberg & Cole 
1992). 

It is hoped that the three methods described here span 
a broad range of plausible biasing prescriptions that might 
occur in the real universe, and that, if we can develop tech- 
niques that are unaifected by which of these schemes is cho- 
sen, these techniques will be insensitive to the way in which 
bias comes about. Further details about biasing schemes will 
follow in Cole et al. (in preparation). 

The density of the biased catalogues is chosen to be 
one galaxy per 25(/i~'^ Mpc)^, roughly four times 0», the 
observed number density of i* galaxies (eg. Peebles 1993, 
§5;Ratcliffo ot al. 1997b). This high sampling density has 
been adopted as we are interested in investigating system- 
atic effects on the redshift-space distortions, and so we want 
to minimize the statistical uncertainties produced by sparse 
sampling. 



3.2 Zel'dovich approximation simulations 

We attempt to model the redshift-space distortions in our 

A'^-body simulations using a Monte Carlo implementation of 
the Zel'dovich approximation. Our Monte Carlo realizations 
have the same 192^ particle grid and initial density field pa- 
rameters as the corresponding A'^-body simulations. Exactly 
the same method is used to perturb the particles as that 
employed to got the initial displacements in the A^-body 
simulation, only with much larger perturbation amplitude 
corresponding to z — rather than the starting redshift 
of the simulation. For the A'^-body simulations we populate 
all the modes in the 192^ fc-space grid including those in 
the corners of the cube. However for these realizations of 
the ZA, we make our treatment consistent with the analytic 
calculation that we outline in Section 4.3, by smoothing the 
density field in the same way, ie. by truncating the power 
spectrum isotropically at /cNq ~ 1.743/). Mpc"'^. 

These realizations of the ZA can be biased using exactly 
the same prescriptions that we use for the A''-body simula- 
tions. The catalogues thus created are then subjected to the 
same multipole analysis as our AT-body simulations. 



3.3 Estimation of multipole moments 

We ensure that the distant observer approximation is satis- 
fied by assigning redshift-space positions to galaxies based 
on their distance along the x-axis of the simulation. Effec- 
tively the simulation box is placed at infinity and the line 
of sight aligned with the a;-ax;is. To minimise the noise in 
our estimates we repeat the analysis with the simulation ro- 



tated to align the line of sight with the y- and then z-axes 
and average the results. 

The galaxies in the simulation are assigned to a 192^ 
grid using the cloud-in-cell (CIC) weighting scheme (Efs- 
tathiou et al. 1985). The density field thus produced then 
has the average galaxy density subtracted from it and is di- 
vided by this average density. The resultant grid is thus the 
fractional overdensity field, 6{r), and the grid is fast Fourier 
transformed to got (5(k). The effect of the CIC assignment is 
to doubly convolve the density field with the shape of a grid 
cell. Thus to deconvolve we divide by the square of Fourier 
transform of the CIC window function. The power spectrum 
estimator is then (5(k)(5(k)*. 

For the biased galaxy simulations a constant equal to 
the inverse of the number density of galaxies in the catalogue 
is subtracted at this point to account for the shot noise in- 
troduced by the Poisson sampling of the galaxies. This is 
not necessary in the unbiased cases as here we use all the 
particles rather than a Poisson sample. 

We now divide fc-space in spherical shells of width 
Afe = 27r/I/box. For each shell we estimate the multipole mo- 
ments of the power spectrum by approximating the integral 
of fi in (2.7) by a discrete sum over modes. This method 
works perfectly well at high-fc, but for the first few shells 
where the sampling of fc-space results in poor sampling in /x 
it leads to a small systematic error in P2{k)/Po{k). An al- 
ternative method of estimating P2{k) and Po{k), which can 
readily be employed for low-fc, is to make a least squares fit 
to the /n dependence of P(fc, n) using a basis of the Legendre 
polynomials Lo,L2 and L4. This second method produces a 
less noisy estimator but still suffers from a small systematic 
error due to the poor sampling in fi. We therefore apply a 
small (5% for the lowest wavenumber), empirically derived 
correction to our first five bins in fc to compensate for this. 

The lowest fe-mode probed using this method is fe = 
(27r/_Lbox), ic. when one wavelength fills the whole box. The 
highest is the Nyquist frequency, fcNq = (27r/Lbox)(Af/2), 
where A^ is the number of cells along one side of the simu- 
lation box. We do not expect our treatment to work all the 
way up to the Nyquist frequency, since the deconvolution 
correction can never recover all the information lost by as- 
signing particles to a discrete grid. This effect can be seen 
in Fig. 1 where the power spectrum turns up for the very 
highest fc-values. For this reason we limit our analysis rather 
conservatively to modes with k < ifcNq- 



4 ANALYTIC MODELS 

If it is assumed that the velocity field, v(r), is generated via 
gravitational instability, then it can be related to the un- 
derlying mass density field. Here we consider three analytic 
models for the relationship of the velocity field to the density 
field, which lead to quantitative descriptions of the redshift- 
space distortion of galaxy clustering. The first is base purely 
on linear theory. In the second, the effect of non-linear ve- 
locities in galaxy clusters is modelled by adding a random 
velocity dispersion to the velocity predicted by linear the- 
ory. The magnitude of this dispersion, CTi, , is not predicted 
by the model, but is expected to bear some relation to the 
typical thermal velocity of galaxies within groups and clus- 
ters. The third model is the Zel'dovich approximation (ZA) 
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Figure 2. Model predictions for the redshift space power spectrum quadrupole to monopole ratio P2{k)/Po{k). The two sets of curves 
are for two unbiased (b = 1) CDM models with T = 0.25 power spectra. The upper curves are for Qq = 1-0, as = 0.55, while the lower 

model is Qq = 0.3, Aq = 0.7 and o-g = 1.13. In both cases the curves show the constant ratio predicted by linear theory (dotted), the 
analytic Zel'dovich Approximation result (solid), and examples of the dispersion model (dashed) with ay = 300km s^^ for Qo = 1.0 
and cTy = 500km s^^ for Qq = 0.3. The error bars show the result of averaging forty random Monte Carlo realizations of the Zel'dovich 
Approximation. These error bars indicate the standard deviation of the sample, and thus are typical of the error expected in a single 
realization. 



(Zel'dovich 1970). This model provides a complete descrip- 
tion of the evolved density and velocity fields, which is ex- 
pected to be accurate in the quasi-linear regime where the 
overdensity is of order unity. Thus this model has the poten- 
tial to give an accurate physical description of redshift-space 
clustering. These three models are described in detail below 
and compared in Fig. 2. In each case we assume that the 
distant observer approximation is well satisfied and that we 
are dealing with volume limited galaxy samples. The eff'ects 
of relaxing these assumptions are discussed in Kaiser (1987). 

4.1 Linear Theory 

In the regime of linear theory Kaiser (1987) showed that the 
redshift-space power spectrum is related to the real space 
power spectrum in a very simple way 



P(fe,/x) = P(fe)(l + /3M' 



13 = f{no)/b : 



(4.1) 



The function /(fio) ~ f^o^ is the logarithmic derivative of 
the fluctuation growth rate (see Peebles 1993, §14; Bouchet 
et al. 1995). The bias factor, h, is an assumed constant relat- 
ing fluctuations in the gala^xy density to those in the mass. 
Applying (2.7) for the I = and I = 2 modes results in 

^ ^ 4/3/3 + 4/3V7 

Po 1 2/3/3 + /3V5' ' 
Thus linear theory predicts a constant quadrupole-to- 
monopole ratio, P2/P0, independent of scale. This model 
is shown by the dotted lines in Fig. 2. 



Here we highlight some of the steps in this derivation 
so that we can compare its assumptions with those of two 
other proposed models. In linear theory the velocity and 
gala^xy density fields are related by the first order continuity 
equation, 



V.v(r) = -f36{r). 



(4.3) 



Considering the density field as a sum over Fourier compo- 
nents, 

5(r) = ^5(k)e-^''-, (4.4) 
the resulting linear theory prediction for the velocity field is 
v(r) = -/3^g5(k)e-*-, (4.5) 



where we have assumed a linear bias between the galaxy and 

mass density perturbations, S{r) = 6(5niass (r) • This leads to 
the following expression for the line-of-sight velocity gradi- 
ent. 



— = -/3p 5(r) 



(4.6) 



The results (4.1) and (4.2) are valid only so long as each of 
the four following constraints is satisfied, so that the corre- 
sponding form of non-linearity can be ignored (Cole et al. 
1994, section 2.2); 

1. dispersion: The velocity dispersion in virialized sys- 
tems is sufficiently small that it can be ignored, i.e. fccr^ <C 1 
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2. dynamical: The linear relationship between velocity 
and overdcnsity (4.5) is accurate. 

3. gradient: Second order terms in the gradient of the 
line-of-sight velocity field can be ignored, i.e. dU jdr <C 1 

4. contrast: Second order terms in the galaxy overden- 
sity can bo ignored, <IC f 

All these constraints will be satisfied on the very largest 
scales, but, depending on the present-day amplitude of 
galaxy clustering and on the values of f2o and fe, we expect 
some or even all of them to be violated on scales of less than 
lOOft-^ Mpc. 

4.2 Linear Theory Plus Dispersion 

The velocity dispersion of galaxies in galaxy clusters is 
typically 800km s~^ (eg. White et al. 1993). Thus, from 
the the constraint that kov ^ 1, we would expect the 
linear theory result (4.1) to apply only on scales where 
fc <C 0.125ft Mpc^^. Since wo arc interested in clustering 
on scales from 10 - 200/i"^ Mpc (0.03 < fe < 0.6/iMpc"^), 
an accurate model will have to take into account these non- 
linear velocities. 

On small scales the redshift space correlation function 
has boon found to be well modelled as a convolution of the 
real-space isotropic correlation function with an exponential 
probability distribution function for line-of-sight velocities 
(Bean et al. 1983; Davis & Peebles 1983; Fisher et al. 1994). 

Park et al. (1994) point out that, since the convolution 
in r-spacc is equivalent to a multiplication in fc-spacc, the 
power spectrum is multiplied by the square of the Fourier 
transform of the velocity distribution function. Peacock & 
Dodds (1994) show that the effects of linear clustering and 
this model of small-scale velocity dispersion can be combined 
to give a redshift space power spectrum 

P(fe,/.) = P(fc)(l + /3MY(l + ^^^) (4.7) 

For this expression it is considerably harder to perform the 
integration over /x required in (2.7) than for the linear theory 
case, but an analytic expression for P2{k) / Po{k) can easily 
be obtained with the aid of mathematics packages such as 
Maple. The rather lengthy formulae are given in Cole et al. 
(1995) section 2.1. 

This model extends linear theory by relaxing the dis- 
persion constraint listed in the previous section, but still 
implicitly assumes the other three. The model of the small- 
scale velocity dispersion is also simplistic in that it takes no 
account of the fact that the velocity dispersion is correlated 
with the density field, ie. the dispersion is higher in high 
density regions such as galaxy clusters. This means that the 
value of av used in this model is only an effective veloc- 
ity dispersion which depends on how galaxies populate the 
clusters, and thus on the bias parameter, b. Examples of this 
model are shown by the short dashed curves in Fig. 2. 

4.3 The Zel'dovich Approximation 

The linear theory approach to modelling the growth of den- 
sity perturbations is only valid if (5(r) <C 1. The refinement 
to linear theory proposed by Zel'dovich (1970) was to formu- 
late a Lagrangian approach. Here each particle is displaced 



from its original position along a straight line defined by 
the direction of the initial velocity field. In comoving co- 
ordinates r, the final position, is related to q, the initial 
position, by 

r = q + d(q,t), d(q,t) = ^. (4.8) 

The Zel'dovich Approximation (ZA) is expected to 
break down at the stage when shell-crossing occurs {S{r) ~ 
1). This is because, under this model, the particles pass right 
through caustics as they continue to move in the direction 
of their original velocity. In contrast, in iV-body simula- 
tions non-linear off'ccts cause the particles to behave as if 
they were 'sticky', and galaxies congregate in high density 
shells or walls. If the power spectrum is not truncated or 
filtered at high spatial frequencies, shell-crossing occurs on 
small scales, with the result that small-scale and some large- 
scale power is erased, and the degree of anisotropy is lower 
than expected. We choose to smooth the initial density field 
for the Zel'dovich realizations on small scales by applying a 
sharp cut-off to the power spectrum at kr = 1.743/i Mpc~^, 
which corresponds to the Nyquist frequency of the grid used 
in the TV-body simulations. 

Fisher & Nusser (1996, henceforth FN96) and Taylor & 
Hamilton (1996) obtained an analytic result for the redshift- 
space distortion produced with the ZA. Following their work 
we have developed an integration technique using an Euler 
method (Press et al. 1992) to deal with diverging oscillations 
of the integrand. Over the range of k we are interested in, the 
method is quite stable, and avoids the complexity described 
in the appendix of Taylor & Hamilton (1996). However, our 
technique breaks down at a scale k ^ O.G/irg. At this point 
oscillations in the radial integrand become too rapid and of 
too great an amplitude for the method to cope. 

This analytic calculation of the quadrupole ratio 
P'zik) / Paik) is compared with the results of averaging many 
Monte Carlo realizations of the ZA in Fig. 2. The close 
agreement of the two methods is very reassuring. It not 
only demonstrates the accuracy of our numerical integration 
technique, but also our implementation of realizations of the 
ZA and the entire procedure of estimating P2(k)/Po{k) from 
particle distributions. 

In Fig. 2 the P2/P0 ratio of the fio = 0.3 model reaches 
a minimum at fe « 0.57iMpc~^ and then slowly increases. 
The ratio behaves in a similar manner for fio = 1.0, but 
with a minimum at k ~ 0.95ft Mpc^^, which lies ofi' the 
right-hand side of Fig. 2. We attribute this behaviour to the 
breakdown of the ZA at scales on which substantial shell 
crossing has occurred. This occurs on larger scales in the 
Q.0 — 0.3 model as it has the larger amplitude of density 
fiuctuations with as ~ 1.13 compared with ag = 0.55 in the 
f2o = 1-0 model. In reality one expects the ratio P2{k) / Po{k) 
to become increasingly negative on small scales due to the 
fingers-of-god produced by the random velocities of galaxies 
in virialized groups and clusters. 

In our implementation of the ZA we have applied very 
little smoothing of the initial density field. We have sim- 
ply truncated the input power spectrum at the Nyquist fre- 
quency of our particle grid used in the Monte Carlo realiza- 
tions described in Section 3.2 to facilitate inter-comparison. 
However, Melott, Pellman & Shandarin (1994) have shown 
that, at least for some statistics, the correspondence between 
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the ZA and N-body simulations can be improved by more 
severe smoothing of the initial density field. They define a 
non-linear wavenumber fcni by 

4tt / P{k)k^di = l (4.9) 
Jo 

and then smooth the initial density field with a Gaussian 
window, cxp (— fc^/2fcs ), with fcs ~ pk-ni- Tlicy find that the 
best choice of the parameter p depends only weakly on the 
shape of the power spectrum and is p ~ 1.25. When com- 
paring to the A'"-body simulations in Section 5 wc investigate 
whether this procedure improves the redshift space predic- 
tions of the ZA. 

The ZA has several advantages over pure linear theory. 
Firstly, equation (4.5) is modified such that the velocity v(r) 
is given in terms of the density field at the initial position q, 
ie. r is replaced by q in (4.4), and on the right-hand-side of 
(4.5). This dynamical relation remains quite accurate until 
shell crossing occurs, whereas the linear theory relation is 
only accurate for 5{r) -C 1. In contrast to the linear theory 
derivation, no further assumptions regarding either the am- 
plitude of the velocity gradient nor the density fiuctuations 
are made in deriving the resulting power spectrum. Thus in 
this respect it is able to deal with both gradient and con- 
trast non-linearity. For modelling redshift-space distortions 
the main shortcoming of the ZA is that it does not model the 
random velocities produced in non-linear relaxed structures 
- instead it produces its own velocity dispersion on small 
scales due to the shell crossing. Thus in this respect it does 
not deal with the dispersion non-linearity. Also in the form 
described above it docs not explicitly deal with the cfi^ccts 
of biasing. Fisher & Nusser (1996) claim that to a good ap- 
proximation bias can be modelled by simply replacing /(fio) 
with (3 = f{i}o)/b - we investigate this in Section 5. 



4.4 Fitting the models 

Our objective is to illustrate the systematic difiierences be- 
tween the A'^-body results and the models rather than to 
propose a methodology to estimate model parameters from 
real estimates of the redshift-space distortion. In order to an- 
alyze the data from our simulations on quasi-linear scales, 
we will perform fits for the dispersion and ZA models to 
the A''-body points and restrict the range of k to the regime 
over which each model can reasonably match the data. 

The dispersion model provides a functional form for the 
shape of the quadrupole estimator, but says nothing about 
the clustering itself. Thus, in fitting to the data, no assump- 
tions need be made about the shape or amplitude of the 
underlying power spectrum, we simply do a minimizar 
tion using /3 and a„ as fitting parameters. 

In contrast, the ZA is much more physically motivated 
in that it attempts to model both the clustering and the 
associated velocity field from first principles. We first at- 
tempt to utilize a simple fitting formula for the shape of 
P2{k) / Po{k) presented in FN96. This formula depends on /3 
and the zero-crossing of the quadrupole, kn\- We find that 
using this approximation introduces a significant offset be- 
tween the best-fit and true values of (3. We therefore opt to 
evaluate the full expression derived by FN96 for the expected 
quadrupole to monopole ratio at each stage in making the 



minimum-x^ fit. This expression gives the distortion as a 
function of /3 and the power spectrum shape and normal- 
ization (r and as). Here we simply adopt the values of F 
and as used in the simulations and use /3 as the single free 
parameter in the model. For real observations it is hoped 
that one could simultaneously fit for these parameters using 
estimates of the redshift-space power spectrum, Po{k). 



5 RESULTS 

5.1 Unbiased Models 

In Fig. 3 we compare the model predictions for the 
quadrupole ratio, P2{k)/ Po{k), with the results of two unbi- 
ased A'^-body simulations. The points with error bars show 
the AT-body results. The model that produces the lower 
quadrupole ratios has flo = 0.3, Ao = 0.7 and as = 1.13, 
and that with the stronger quadrupole signal has fio = 1 and 
as = 0.55. The error bars placed on the iV-body results are 
indicative of the statistical error we expect on an estimate 
of the quadrupole ratio, from a single realization, and are 
obtained by taking the standard deviation of forty Monte 
Carlo realizations of the ZA. It is clear that, apart from the 
first few points with k < 0.07/iMpc~^, the deviation from 
the constant ratio predicted by linear theory is large. 

The solid lines in Fig. 3 are each for one ZA realization 
starting from the same initial density field as the correspond- 
ing A-body simulation. This model reproduces the A-body 
results extremely accurately at low-fe. For the Qo = 0.3 
case the ZA and iV-body results match well up to fe « 
0.25/iMpc~^ which is beyond the scale at which P2{k)/Po{k) 
becomes negative. This result is in complete agreement with 
the findings of FN96. The only A-body results they studied 
were for a similarly normalized Qo = 0.3,Ao = 0.7 model. 
In the case of ilo = 1 we find much less impressive agree- 
ment between the ZA and A-body simulations. The results 
begin to diverge at fc > 0.15/iMpc~^ and differ very signifi- 
cantly at A; = 0.425/1 Mpc"\ where P2ik)/Po{k) = for the 
A-body simulation. 

In Fig. 4 wc show another comparison between these 
A-body and ZA realizations. Here we show the ratio of esti- 
mates of their final power spectra to the initial power spectra 
evolved to z = assuming linear theory. We see that in the 
A-body simulations the non-linear growth of structure leads 
to an increase in power on small scales. This is strongest 
for Oo = 0.3 which has the highest as and is therefore the 
most non-linear. In contrast the behaviour of the ZA simu- 
lations is the exact opposite - shell crossing erases power on 
small scales. In this respect the agreement between ZA and 
A-body is particularly poor for the f2o = 0.3 model. Even 
at mildly non-linear scales this discrepancy is quite severe, 
which suggests to us that the success of the ZA at modelling 
the ratio P2{k) / Po{k) should be treated with skepticism. 

The possibility of improving the correspondence with 
the A-body results by smoothing the initial density field 
before applying the ZA, as in the truncated ZA of Melott 
et al. (1994), is investigated in Fig. 5. The introduction of 
additional smoothing worsens the agreement between the A- 
body results and those of the ZA. In particular for the degree 
of smoothing advocated by Melott et al., corresponding to 
p m 1.25, the agreement between AT-body and ZA is very 
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Figure 3. The quadrupole ratio P2{k)/Po{k) for two unbiased AT-body simulations compared with the corresponding Zel'dovich approx- 
imation simulations and two model fits. The upper curves are for Qq = 1 and as = 0.55 and the lower curves for Qq = 0.3, Ao = 0.7 
and erg = 1.13. The error bars on the Af-body results arc indicative of the statistical error in our estimates and are estimated from the 
standard deviation between forty ZA realizations. The solid lines show the results from single ZA realizations starting from the same 
density fields as the A^-body simulations. The long dashed and short dashed curves are model fits for the ZA and dispersion model 
respectively. Only the first ten points have been used in constraining the ZA fits, and the first thirty for constraining the dispersion 
model. The corresponding best fit parameters and confidence intervals are shown Fig. 6. 
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Figure 4. The ratio of estimates of the real-space evolved power spectra, P{k), to the linear theory predictions, PLlNEAR(fe), for the 
Af-body simulations and the corresponding ZA simulations. The models are the same f2o = 1 and f2o = 0.3 models studied in Fig. 3. 
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Figure 5. The dependence of the Zel'dovich approximation prediction for the ratio P2{k)/Po{k) on the extent to which the initial 
density field is smoothed before applying the ZA. The error bars show the A^-body result for the Oo = 1-0 model taken from Fig. 3. The 

dotted curve shows the ZA result for the minimal amount of smoothing, ic, truncation at the Nyquist frequency. The remaining curves 
show the effect of gradually increasing the smoothing using an additional Gaussian smoothing with ks = pk^\ (see equation 4.9) and 
p = 5.0,2.0,1.25 and 1.0. 
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Figure 6. The 68% and 99.7% confidence intervals (1- & 3-ct) of the model parameters for the fits shown in Fig. 3 and Fig. 7. The ellipses 
are for the two-parameter dispersion model for the three different simulations, and the horizontal lines for single-parameter ZA model 
for the two unbiased models for which this provides a reasonably good fit. The best fit values of /? and the per degree of freedom for 
the best fit models are given in Table 1. 
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poor. Presumably the smoothing being apphed is washing 
out the distortions we are trying to measure. However, we 
find that in real space a mild amount of smoothing can pro- 
duce a power spectrum that is slightly closer to the iV-body 
result, though this is a rather small effect. 

The ZA and dispersion model fits are shown by the 
smooth long and short dashed curves in Fig. 3. The corre- 
sponding best fit model parameters and confidence inter- 
vals are shown in Fig. 6 and Table 1. Here the horizontal 
lines bracket the 68 per cent and 99.7 per cent (1- and 3- 
cr) confidence intervals of /3 for the ZA fit. The ellipses show 
the corresponding confidence intervals for the two-parameter 
dispersion model. Both ZA fits give a best-fit that is within 
1 — a of the theoretical value, so there seems to be no appre- 
ciable systematic error in our estimate. However, this result 
is achieved by limiting the range of the fit to the first ten 
bins in fc-space, ie. k < 0.182/i Mpc~^, in both models. At 
smaller scales than this, particularly in the flo = 1 model, 
the fits can be seen to strongly overestimate the redshift- 
space distortion. Thus if a fit were made to data covering a 
wider range of scales, /3 would be severely underestimated. 
In contrast, the dispersion model fits remain in reasonable 
agreement with the data up to much larger k. We fit our 
data up to the thirtieth data bin, corresponding to a length 
scale of 11.5/i~^ Mpc. For flo = 0.3, the dispersion model 
analysis seems to be offset towards a low estimate of f5, but 
this effect is fairly small given the size of the error contours. 
The fit itself is a very good over the whole range of k in 
question. However, for Qq = 1, (3 is systematically under- 
estimated by more than 10 per cent at a high level of sig- 
nificance. This offset can be understood by reference to the 
models shown in Fig. 2. In comparison to the ZA, which is 
an accurate fit at low-fc, the dispersion model curve changes 
slope too rapidly and is thus too shallow at low-/c. Thus, 
when data at low-fc are used to constrain the fit, the value 
of /3 is underestimated. We found this behaviour to be the 
same for a variant of the dispersion model in which the small 
scale random velocity is assumed to be Gaussian rather than 
exponentially distributed. This alternative Gaussian disper- 
sion model produces a worse fit overall, and a more biased 
estimate of l3. The same effect is also expected to introduce 
an offset in the f^o = 0.3 simulation, but it is rather smaller 
for this case and so is not easily detected. 

As shown in Table 1, we have calculated the 1- 
dimensional rms peculiar velocity of particles in each simu- 
lation. af^\ in an attempt to look for a physical significance 
of the av value found in the dispersion model fit. It is clear 
from these measurements that there is no simple relation- 
ship between the two quantities over the range covered by 
our simulations. 

5.2 Biased Models 

In Fig. 7 we show the model predictions for the quadrupole 
ratio, P-zlk) / Po{k) , for a biased model with Slo = 1-0. The 
heavy curves show the results of selecting biased galaxy cat- 
alogues by the three methods described in Section 3.1. The 
differences between the models at large scales are purely sta- 
tistical, since different galaxies have been selected in the dif- 
ferent samples. At high fc the two Lagrangian methods agree 
well, but the Eulerian one systematically shows less distor- 
tion. This is to be expected as this method populates evenly 
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Table 1. Best fit /3 and goodness of fit (x^ per degree of freedom) 
indicator for our simulations under the Zel'dovich and dispersion 
models. For the dispersion model we also show the best fit values 
of (T„ and compare with the 1-dimcnsional velocity dispersion of 
the simulations, c^"^. Error contours for the fits are shown in 
Fig. 6. 

all regions with density higher than the threshold, rather 
than giving extra weight to the areas with higher density 
(and thus velocity dispersion), as the other two models do. 

The long dashed curve shows the prediction for the an- 
alytic ZA where, as suggested by FN96, the model adopted 
assumes that galaxies are unbiased, u% = crf^', but with 
no** = 1/1.81 to match the value of /3 = 1/6. This curve 
should be compared with the results from a biased catalogue 
constructed from a realization of the ZA which is shown by 
the thin solid curve. Clearly these two curves do not agree. 
Thus in the ZA, bias cannot be taken account of by merely 
modifying the value of /(fio) as suggested by FN96. This 
is perhaps to be expected as bias will not only boost the 
galaxy density fluctuations with respect to the underlying 
mass distribution, but will also preferentially place galax- 
ies in dense structures where random velocities produced by 
virialized structures are largest. The small scale velocity field 
produced by shell crossing in the ZA is not a good model of 
the random velocities produced in non-linear regions. It is 
for this reason that the biased ZA realization is a poor fit 
to the A'"-body results. It diverges from the A'"-body result 
more rapidly than the corresponding model for the unbiased 
fio = 1 simulation. 

The short dashed curve is a fit of the dispersion model to 
the iV-body peaks biased galaxy catalogue. This model, with 
an exponential velocity distribution, determines (3 very ac- 
curately and is a good fit to the data over this range of scale. 
The fitted velocity dispersion is ay = 380km s~^, 70km s~^ 
higher than the unbiased version of the same iV-body simula- 
tion, which illustrates how bias preferentially places galaxies 
in environments with higher thermal velocities. 



6 DISCUSSION 

We have used a set of high resolution A/^-body simulations to 
investigate the accuracy of two models of the redshift-space 
distortion of galaxy clustering on scales of 10-200/i~^ Mpc. 
We conclude that neither the Zel'dovich Approximation 
(Fisher & Nusser 1996; Taylor & Hamilton 1996) nor the lin- 
ear theory plus random velocity dispersion model (Peacock 
& Dodds 1994) provide an accurate description of redshift- 
space clustering, lu general both models could lead to signif- 
icant systematic errors in the estimation of /J = /b when 
appUed to the high precision data that will be available from 
the large 2dF and SDSS galaxy redshift surveys. 

The Zel'dovich Approximation provides a good fit in the 
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Figure 7. The quadrupole ratio, P2{k)/Po{k), for biased fJo = 1 galaxy catalogues with b = <t|* /as = 1.81, /3 = 0.55. The AT-body 
lines show the results for our three different methods of constructing biased galaxy catalogues from the iV-body simulation, as outlined 
in Section 3.1. The thin solid line shows the result of a ZA realization that has been biased in the same as one of the A^-body catalogues. 
The long dashed curve shows the prediction of the ZA made by assuming eg = o-f and = /3 as proposed by FN96. The short 
dashed line is a fit of the dispersion model to one of the AT-body catalogues. 



unbiased cases on the very largest scales that we have inves- 
tigated. However, only for the low Qq model that we exam- 
ined does this model remain accurate over the full range of 
scales for which we mesisure a positive quadrupole distortion 
(P2{k) / Po{k) > 0). In the case of biased galaxy catalogues 
the ZA is a very poor description of redshift clustering apart 
from on the very largest scales. Smoothing the initial den- 
sity field prior to applying the ZA, as in the Truncated ZA 
(Melott, Pellman & Shandarin 1994), weis found to further 
worsen the agreement with the iV-body results of redshift- 
spacc distortions. In general, then, we would only advocate 
using the ZA for analysis of unbiased universes on scales 
down to ~ 40/i"^ Mpc. 

The two-parameter dispersion model was found to pro- 
duce more acceptable fits to the iV-body results. However, 
only in the cases in which the small scale velocity dispersion 
was large did these fits yield values of /3 that were not sig- 
nificantly biased with respect to the true values. This model 
also yields a value for the velocity dispersion parameter , 
but this depends on the degree and form of galaxy bieis and 
so is difficult to directly relate to an interesting physical 
quantity. It should be noted that, although the dispersion 
model was unable to deal with the unbiased Qo = 1.0 simu- 
lation, and in fact extracted a value of (3 offset from the true 
value by several-sigma, results from a variety of cosmolog- 
ical methods arc making it seem increasingly unlikely that 
we live in such a universe - constraints from the observed 
as and current /3 estimates require a low-Oo or biased cos- 
mology. Thus the dispersion model may well prove to be a 



useful tool for extracting information from redshift surveys 
on intermediate scales. 

There appear to bo two distinct ways to proceed to 
remedy or by-pass the above short-comings. 

1) Improve the analytic models so as to produce ac- 
curate predictions over the full range of scales that will be 
probed by the redshift surveys. 

A promising approach is an extension of the model dis- 
cussed in Fisher (1995), where linear distortion effects on 
^(fj, tt) are modelled as a convolution with a Gaussian ve- 
locity distribution with scale- and orientation-dependent ve- 
locity dispersion. Fisher suggests that this model may be 
adapted to deal with distortions in the highly non-linear 
regime by replacing the Gaussian with a generalized distri- 
bution function which is Gaussian on large scales but tends 
to an exponential with isotropic and constant velocity dis- 
persion on small scales. 

Hamilton et al. (1991) developed an accurate method of 
obtaining the non-linear correlation function in real space, 
^(r), using a universal scaling relation. Mo, Jing, & Borner 
(1996) outline a technique by which the quantities which 
define the velocity field (ie. mean pairwise peculiar velocity, 
Vi2{r), pairwise pecuhar velocity dispersion, (ui2 (?■)), and 
mean square peculiar velocity, {vf)), can be modelled for a 
given cosmology and initial power spectrum. 

Given reliable models for these quantities. Fisher's 
treatment should provide a way of translating ^(r) into the 
redshift space statistic ^{a, n), and thus calculating precisely 
the anisotropy expected on any given scale. However, this 
method is still subject to problems - namely, finding out 
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Figure 8. A compajrison of the quadrupole ratio, P2{k) / Po{k) before and after identifying and collapsing the galaxy clusters. The dashed 
lines show the original shape for the three simulations, where we have used the peaks biasing method for the biased case. The solid lines 
the result after cluster collapse for the three models. We also plot the other two variants of the biasing model in the collapsed case, using 
a dotted line for method 2 and dot-dashed for method 3. 



exactly how the velocity distribution function behaves on 
intermediate scales, and developing a way for modelling the 
effects of galaxy bias. We intend to investigate this in a fu- 
ture paper. 

2) As our computers increase in memory and CPU 
power, the 'brute force' approach of simply running iV-body 
simulations for the whole range of parameter space becomes 
increasingly plausible. Here the dispersion model may be 
very useful in fitting and suggesting the approximate value 
of f3, and thus putting restrictions on the combinations of 
parameter values worth investigating. 

If one is resorting to simulations then it may be possi- 
ble to filter the data to remove some non-lineaxities before 
commencing the analysis of redshift-space distortions. An 
interesting procedure for doing this is to identify and col- 
lapse galaxy clusters. Since their high velocity dispersion 
is responsible for much of the suppression of P2/P0, their 
removal might lead to more robust estimates of /3. Fig. 8 il- 
lustrates the effect of collapsing clusters in our simulations. 
As expected, the ariisotropy is much closer to the linear the- 
ory prediction on large and intermediate scales. On small 
scales the quadrupole ratio converges on zero rather than 
becoming strongly negative. Cluster membership was estab- 
lished using a friends-of-friends algorithm, in which any pair 
of galaxies with separation less than 0.2 times the mean sep- 
aration are classed as being in the same cluster. Each cluster 
member then has its velocity replaced with the mean veloc- 
ity of the cluster itself, effectively eliminating the internal 
velocity dispersion of these virialized structures. 

The identification of galaxy clusters will be much more 
complex when dealing with real galaxy redshift surveys. 



Moore, Frenk & White (1993) describe a way of apply- 
ing this technique to the CfA survey using a redshift- 
and orientation-dependent linking length to define groups. 
Whatever algorithm one adopts can always be tested by run- 
ning it on mock galaxy catalogues made from the iV-body 
simulations. 

Despite the failings of existing models pointed out in 
this work, redshift-space distortion analysis remains very 

interesting and promises to be of great use in extracting 
information from the next generation of galaxy surveys. 
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